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We calculate the fermion self-energy at 0{as) for the Symanzik improved staggered fermion action 
used in recent lattice simulations by the MILC collaboration. We demonstrate that the algebraic 
approach to lattice perturbation theory, suggested by us recently Q, is a powerful tool also for 
improved lattice actions. 
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I. INTRODUCTION 

In a recent paper P] , we have shown that methods fa- 
miliar from continuum perturbation theory can be used 
to perform calculations in lattice perturbation theory an- 
alytically, except for the numerical evaluation of a small 
number of lattice "master" -integrals. In short, the cal- 
culation proceeds as follows: in a first step, a given loop 
integral is expanded around the continuum limit using 
the technique of "asymptotic expansions" for Feynman 
diagrams, familiar 2] from continuum perturbation the- 
ory. The expansion splits the original integral into con- 
tinuum integrals and a set of lattice tadpole integrals. 
These tadpole integrals are independent of external mo- 
menta and masses and are interrelated through a set of 
algebraic identities. Using these relations, the tadpole in- 
tegrals can be expressed through a small number of mas- 
ter integrals which are evaluated numerically. In jj, we 
have illustrated this strategy by evaluating the gluon and 
fermion one-loop self-energies with the standard Wilson 
action. 

We have also asserted that the above techniques are 
not limited to a particular lattice action or to one loop. 
While correct in principle, in practice this could turn 
out to be of little use, given the limitations of available 
computing power. The goal of the present paper is to 
investigate if the method proposed in remains useful 
for improved lattice actions which are employed in many 
lattice simulations. In what follows, we perform an an- 
alytic one loop calculation with the lattice QCD action 
used by the MILC collaboration in their recent simula- 
tions with three dynamical quarksH • This is a one- loop 
Symanzik improved gluon action |j] coupled to improved 
Kogut-Susskind fermions This action is accurate 

up to errors of 0{a?'as,a^). We find that the amount 
of computer algebra needed for perturbative calculations 
with such an action is significantly larger than for the 
simple Wilson action, but they nevertheless turn out to 
be feasible. 

Let us stress an important difference between the 



purely numerical approach used e.g. in and the al- 
gebraic method Q which we use in this paper. In order 
to make use of all algebraic relations, not only the inte- 
grals that appear explicitly in the expansion of a given 
graph, but the entire class of lattice tadpole integrals rel- 
evant for a given lattice action needs to be considered. 
For the present calculation we had to produce a database 
of roughly SO'OOO tadpole integrals that are all expressed 
in terms of 16 master integrals. This took about a week 
of CPU time running Maple on a desktop computer. Less 
than 2000 of those integrals actually occurred in our dia- 
grams. On the other hand, once the database is set up, it 
can be used for any one-loop calculation involving stag- 
gered fermions. In this paper we deal with the simplest 
example of the quark self-energy and derive the relation 
between the bare and the pole mass of the quark at or- 
der 0{as). We plan to return to other applications in 
the near future. 

This paper is organized as follows. We start by sum- 
marizing the method proposed in and by pointing 
out some peculiarities in its application to staggered 
fermions. The presence of the fermion doublers mani- 
fests itself in the appearance of multiple "soft" regions. 
We then explain our treatment of the improvement terms 
in the action. In our calculation we expand the im- 
proved gluon and fermion propagators in a geometric 
series around the standard propagators. Since loop in- 
tegrals are analytic in the improvement terms, this ex- 
pansion can be performed at the level of the integrand. 
As a matter of principle, it can be carried to any level of 
precision desired. Finally, we present the result for the 
fermion self-energy and conclude. 



II. ASYMPTOTIC EXPANSION AROUND THE 
CONTINUUM LIMIT 

In this Section we briefly describe the main idea of our 
approach to lattice perturbation theory 0. We are in- 
terested in the calculation of quantities in the continuum 
limit, where the particle masses and external momenta 
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are much smaller than the inverse lattice spacing. It is 
possible to construct the expansion of a lattice integral 
around the continuum limit by expanding the integrand, 
at the expense of introducing an additional (analytic) 
regulator. 

It turns out to be convenient to map the integra- 
tion region — tt < fc^ < tt to infinite volume, by changing 
the integration variables to 

7?^ =tan(A:^/2). (1) 

Furthermore, we introduce an intermediate regulator by 
adding a small quantity S to the power of one of the prop- 
agators. Keeping the regularization parameter S non- 
zero, we then expand the integrand of the lattice loop 
integral in the so-called soft and hard regions. The orig- 
inal integral is given by the sum of the soft and hard 
contributions; in this sum all the terms that are singu- 
lar in the limit S ^ cancel out. The soft and hard 
contributions are calculated as follows: 

• Soft: assume that all the components of the loop 
momentum 77 are small, i.e. they comparable to 
masses m and external momenta p in the problem, 
r]i m, p <ti I. Perform the Taylor expansion of 
the integrand in the small quantities rji and m,p. 
The expansion coefficients in this region are stan- 
dard continuum one loop integrals, regularized an- 
alytically. No restriction on the integration region 
is introduced. 

• Hard: assume that all the components of the loop 
momentum are large, 77^ ~ 1 ^ m,p and Taylor 
expand the integrand in m,p. The expansion coef- 
ficients are massless lattice tadpole integrals. 

As in dimensional regularization, scaleless integrals are 
set to zero. 

While the soft part is given by continuum integrals, 
the hard part is given by lattice tadpole integrals that 
do not depend on the process under consideration. Intu- 
itively, this is the case because the hard part arises from 
the integration region where the loop momentum is com- 
parable to the inverse lattice spacing and therefore much 
larger than the external momenta. The analytic regu- 
larization which we use in what follows, turns out to be 
very efficient in separating soft and hard momenta modes 
and allows us to treat them separately. It is important to 
realize that since the integrals occurring in the hard part 
are process independent, they can be calculated once and 
for all for a given lattice action. 

In any realistic calculation, many lattice tadpole inte- 
grals appear. Fortunately, not all of them are indepen- 
dent. Relations between the integrals are obtained after 
observing that the integral of a total derivative with re- 
spect to one of the variables r/^ vanishes (there are no 
surface terms, thanks to analytic regularization). By ex- 
plicitly computing the derivative, one finds the algebraic 
relations (recurrence relations) between the lattice tad- 
pole integrals. It is possible to solve these relations using 



the Gaussian elimination method and express all the lat- 
tice tadpole integrals through a small number of master 
integrals. For more details we refer to 

The strategy outlined above does not rely on the spe- 
cific form of the propagator. In 0, we have applied the 
technique to HQET and to QCD with Wilson fermions 
and calculated a number of QCD one-loop self-energies. 
We now show how to use the method to calculate loop in- 
tegrals for staggered fermions. The staggered fermion ac- 
tion is obtained by reducing the number of components of 
the fermion field after spin diagonalizing the naive lattice 
fermion action. In perturbative calculations it is conve- 
nient to work with the naive fermion action and perform 
the staggering only at the end of the calculation. 

The fermion doubling inherent in the naive discretiza- 
tion of the fermion action manifests itself in the appear- 
ance of multiple soft regions. The sixteen zeros of the 
inverse propagator of naive fermions. 



in the Brillouin zone give rise to sixteen propagating 
fermions. Correspondingly loop integrals which involve 
naive fermion propagators have sixteen soft regions. One 
of those is the region where all components of the loop 
momentum r]i are small. The fifteen additional doubler 
contributions arise after transforming one or several com- 
ponents of the integration momentum as rji — > l/rji and 
then expanding the resulting integrand around small 
For purely fermionic integrals, the sixteen soft contribu- 
tions are all equal. In integrals with both boson and 
fermion propagators, this is not the case: since the bo- 
son propagator is far off-shell in all corners of the Bril- 
louin zone, the doubler contributions are different from 
the genuine soft contribution where all the momenta are 
small compared to inverse lattice spacing. 

The contribution of the hard part is obtained as usual, 
by expanding the loop integral in external momenta and 
particle masses. The lattice tadpole integrals that occur 
in the case at hand are 



(lail, n, to) = TT / ^^^5T — (3) 



where Db is the standard bosonic propagator expressed 
through variables "q. Let us stress that the regulator has 
to be on the fermion propagator in order to regulate both 
the singularities at 77^ = and 77^ = 00. The integration- 
by-parts and partial fractioning relations for this class of 
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integrals are 



(4) 



— + 2 Ui (ai — I) + 2 ai {ai — 1) nn 

+2 a^{ai- 1) (2 - 1) (m + 5) mj H 

= a, (1- a,)m\H, 

^ i=i ^ 

= aOnjff. 

i=i ' 



The operator a.^ raises the index by one unit, etc. 
The above relations are sufficient to reduce all integrals 
H to sixteen master integrals. These can be chosen to 
be convergent (in the limit (5 — s- 0) and are calculated 
numerically. An efficient method for their evaluation is 
described in 0. 



III. IMPROVED ACTIONS 



therefore involves 12 powers of the lattice momentum A:^ 
just from the vertex plus additional powers arising from 
the expansion of the two propagators. This is why we 
were forced to compute approximately thirty thousand 
lattice tadpole integrals to have the full reduction of the 
integrals that appear in the calculation of the fermion 
self-energy. 



A. Improved staggered fermions 

The naive discretization of the free fermion action can 
be improved by adding a third lattice derivative term to 
the action: 



The associated fermion propagator 



D/{k) =i^7^sinfc^(l + -sin2 fc^) 



(5) 



(6) 



The lattice artifacts arising in calculations with the 
simplest discretization of the QCD action turn out to be 
uncomfortably large, making it difficult to extract physi- 
cal quantities from simulations with currently accessible 
lattice spacings. For this reason, many simulations are 
done with improved actions, which include higher dimen- 
sional operators whose coefficients are tuned in such a 
way as to reduce the cutoff effects arising at finite lattice 
spacing. 

In principle one could apply the procedure for per- 
turbative calculations described above directly to im- 
proved actions. The higher derivative terms encountered 
in these actions complicate the calculation in two ways: 
(1) The recursion relations for the improved propagators 
are much more complicated than for the standard case 
which might lead to problems with finding their solu- 
tions and, also, to a larger number of master integrals; (2) 
The derivative interactions make it necessary to calculate 
hard integrals with more complicated "tensor" structure, 
i.e. higher powers of a, in Eq.Q. 

In the following we avoid the first complication by ex- 
panding the improved propagators around the standard 
ones. In the soft region, where the momentum flowing 
through the propagators is small, this expansion is built 
into our formalism since for small momenta the higher 
dimensional operators are suppressed by powers of the 
lattice spacing. For large momentum, on the other hand, 
these operators are suppressed only numerically, but not 
parametrically. As we will show this suffices to get the 
result with reasonable precision. 

The disadvantage of expanding the propagators is that 
very high powers of lattice derivatives appear. For ex- 
ample, the quark gluon coupling of the Symanzik im- 
proved staggered quark action involves terms with up 
to six derivatives. The usual fermion self-energy graph 



reproduces, for 6=1, the free continuum propagator 
up to terms suppressed as 0(fc"') at low energies. Note 
that the third order derivative in the above action not 
only improves the fermion propagator at the origin of 
the Brillouin zone, but also in all of its corners where 
some components of the momentum are tt and doublers 
appear. 

From the above expression it is obvious that the im- 
provement term proportional to b is numerically sup- 
pressed over the entire momentum range — tt < fc < tt. 
For the propagator at the origin 
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(2^)4 ^^sinfc2 (i + &sin2fc^)2 
0.6197336 - 0.1267555 b + 0.0214647 b^ 
- 0.00331 b^ + 0.0004803 b^ 4 



(7) 



the expansion in b is roughly a geometric series with an 
expansion parameter of 1/6. 

The above action is no longer free of 0{a?') lattice ar- 
tifacts, once a gauge field is coupled to the fermion. The 
standard discretization of the covariant derivative 
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Ux„ 



(8) 



leads to "flavor" -changing interactions, in which a 
fermion emits gluon with momentum of order tt and scat- 
ters into one of its fifteen partners. The Feynman rule 
for the quark-gluon coupling of naive fermions 

I g {2^f S\k, k, + k,) ^ 7m cos(^^) 
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shows that the flavor-changing and the flavor-neutral in- 
teraction have the same strength. However, in flavor 
changing interactions momentum of the order of the in- 
verse lattice spacing is transferred to the gluon, taking 
it far off-shell and thus suppressing the interaction by 
0(a2). 

To remove the flavor changing interactions at tree level, 
Lepage suggested @ to replace the link field C/^(x) in the 
covariant derivative by a "fat" link 



v 




4 A. 



(9) 

The derivative interactions in the first term of Eq.© 
have been chosen such that the "flavor" -changing inter- 
actions get suppressed, while the second term improves 
the "flavor" -conserving part. Note that the replacement 
U^{x) — > V^{x) is only done for the flrst derivative term 
in Eq. 



B. Improved gluon propagator 

The tree level improved gauge action has the form 

= E {ci^M-'(a^) + C2(i?M^' + ^^'m)}. (10) 

where /3 = 6/(7o, ci — 5/3 and C2 = —1/12. The plaque- 
tte P^i, is the 1x1 Wilson-loop in the /i:/-plane and the 
rectangle i?^^ is the loop with two steps in the /i and one 
in the v direction. Parameterizing the link between the 
lattice site n and n + as 



the quadratic part of the gluon Lagrangian in Fourier 
space is 

Ca2 = A^(-p) [D^, + E^,\ A,{p) (11) 
with [p^ = 2sinp^/2) 



1 r 

12 



1 + e 



(12) 



The first term D^i, corresponds to the usual plaquette 
action obtained by setting ci = 1 C2 = in Ea. (PT1|l . 
The corresponding propagator ispL2| 
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(13) 



We treat i? as a perturbation and expand the improved 
propagator as 

(^imp = G-bG-E-G + b^G-E-G-E-G-... 



The hat variables vary between —2 < p^ < 2. The ef- 
fect of the improvement is largest in the corners of the 
Brillouin zone, where the above expansion is a geomet- 
ric series with expansion parameter 2/3. For the loop 
graphs calculated below, the expansion parameter turns 
out to be around 1/2. To give an example, consider the 
expansion of the /i/i-component of the gluon propagator 
at the origin. It has the following form: 

G]^P{x = 0) = 0.1549334 (1 + e/4) - 0.0366778 b 
+ 0.014347 b^ - 0.0063397 6^ + 0.003045 b* 

- 0.0015511 b^ + 0.0008239 b^ + . . . (14) 

All given digits are significant. To a good approximation, 
the expansion is an alternating geometric series with ex- 
pansion coefficient 1/2. Assuming that the same pattern 
is found at higher orders we estimate the remainder of 
0(6") to be —1/3 of the last computed term and esti- 
mate 



G;,7(0) = 0.1549(1 



^/4) - 0.0266 ± 0.0001. (15) 



The error has been estimated as half of the remainder. 
We will encounter the same behavior for the quark self- 
energy graphs and calculate our final results also in that 
case with the same prescription. 

Let us note that Ea. (|14|l allows us to compute the so- 
called "mean link" , a quantity needed to perform tadpole 
improvement 8] of lattice perturbation theory. At 0{as) 
the mean link in Landau gauge ^ = — 1 is defined as 

{U,{x)/3) = l + ^GlJiO). (16) 

. Using Ea. H14|) we easily find: 

{Uf,{x)/3) = 1 + asUQ = 1 + as( 0.7504 ± 0.001) , (17) 
in agreement with Ref. 0- 

IV. QUARK SELF-ENERGY 

The one loop quark self-energy has the form 

E(p) = as [ij/S^^^ + m E^^)! + o{a^) . (18) 

Below, we separately give the contributions of the two 
graphs in Fig. ^ Since the expressions for the graphs in 
terms of master integrals are neither short nor particu- 
larly illuminating, we have evaluated them numerically. 
The parameter b — 1 indicates the order in the expansion 
of the propagator. 

Soft part 

The leading order soft part is obtained by evaluating 
the continuum diagram in analytic regularization. The 
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FIG. 1: The two graphs contributing to the fermion self- 
energy. 



reason for that is that all the improvements are sublead- 
ing and also the fermion doublers decouple in the con- 
tinuum limit. The doublers contribution is suppressed 
because the appearance of a doubler requires at least one 
of the gluon momenta components to be of the order of 
the lattice spacing. 

We therefore get the same result as for Wilson fermions 



-(1) 

■^soft 

.(2) 
■^soft 



(4 + { 



- + (l - x^) n — 
2 ^ ^ X 



Ls -l + {l + x) In ■ 



1 



with X = w? /p^, Ls = 1/(5 — ln(rn^/4). 



Hard part 

We now compute the hard part of the tadpole and the 
sunset diagrams. As it turns out, the tadpole contribu- 
tion is roughly an order of magnitude larger than that of 
a sunset; for this reason we have to evaluate more terms 
in the expansion for the tadpole than for the sunset. We 
obtain 



.(1) 

■^hard 



.(2) _ 



2.6524- 1.4214 6 + 0.5334 6^ 

-0.2452 6^ + 0.1233 6'^ 

-0.0655 + (0.3245 - 0.0345 b) ^ 

1.60(1) + 0.290^ 

0. 



(19) 



For the hard part of the sunset diagram we find: 



"^hard 



+2) _ 
■^hard ~ 



0.8327 + 0.0880 6- 0.0183 5^ 



3nS 

-(0.6328 - 0.345 6) 



4 + e 
37r5 



- 0.9821 + 0.3279 b - 0.0857 b^. (20) 



The first number is the contribution of the tadpole graph, 
while the second one corresponds to the sunset diagram. 
Note that the hard and the soft contributions to the pole 
mass are not gauge invariant separately. This is the con- 
sequence of putting the regulator S on the fermion prop- 
agator, which we are forced to do to properly deal with 
the doublers. The gauge dependence only cancels in the 
final result for the pole mass, when the hard and the soft 
contributions are added together. 

The effect of the improvement terms in the lattice ac- 
tion is large. This can be seen by comparing with the 
result obtained when all the improvements are switched 
off. In that case the relation between the pole and bare 
masses changes to 



m: 



pole 



= mbarc(l+«. A'r™): 

2 

'-lnm + 4.586. 



(22) 



One sees that the renormalization constant in this case 
is roughly a factor of two larger than for the improved 
action. The reduction arises because the smearing of the 
link in Eq. (|5J| strongly reduces the quark gluon interac- 
tion, while the effect of all other improvements is signifi- 
cantly smaller. 

Note that the result in Ea. (|21|) is not "tadpole- 
improved" 0. For a tadpole improved action, an ad- 
ditional tree level contribution 



A' 



tad 



,(2) 



(23) 



arises, which tends to cancel the tadpole graph. The 

(2) 

quantity uj, is the 0{as) coefficient of the mean link in 
Landau gauge, given in Ea. lfTC|l so that 

^tad ^ _i_688(2). (24) 

Putting everything together, restoring the lattice spacing 
and combining the numerical errors in quadratures, we 
derive the following relation between the pole and the 
bare quark mass in the continuum limit: 



^polc — ^barc 



1 



■ln?7ia + 0.572(14) 



(25) 



Pole mass 



Our results agree with the findings of Ref . i2j . Note that 
this reference uses the mean field uq defined from the 
average plaquette, for which Ug^^ sa 0.7671 7| . 



It is straightforward to extract the pole mass from the 
above formulas. The pole mass is defined as the position 
of the pole of the fermion propagator. It therefore reads: 

Wpolc = = mbarc (1 + as Ai), 

Ai = 

= --lnm+ 1.60(1) + 0.66(1). (21) 



V. CONCLUSION 

Perturbative calculations with lattice regulariza- 
tion are required to provide a bridge between non- 
perturbative simulations and the continuum physics. In 
this paper, we have described an analytic calculation of 
the fermion self-energy of improved staggered fermions 
and the Liischer-Weisz action for gluons. This action has 
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been recently used for non-perturbative lattice simula- 
tions by the MILC collaboration. 

Our main result, the relation between the pole and the 
bare lattice quark masses is given in Ea. (|25l) . We ob- 
serve that the renormalization constant for the improved 
action is rather modest. This is mainly a consequence 
of the smearing of the quark-gluon interaction that was 
introduced to strongly suppress the interaction between 
the fermions in different corners of the Brillouin zone. Let 
us stress that in the current calculation the doublers do 
not contribute in the continuum limit even for the naive 
action; the strong change in the renormalization constant 
therefore has not so much to do with their decoupling as 
with a significant reduction of the quark-gluon interac- 
tion interaction over the entire Brillouin zone. 



This calculation demonstrates that algebraic tech- 
niques for lattice perturbation theory, based on the meth- 
ods developed for continuum perturbation theory, are 
quite efficient. With the large database for one- loop lat- 
tice tadpole integrals for staggered fermions, which we 
have produced in the course of this calculations, other 
one-loop perturbative lattice calculations for improved 
lattice actions seem to be within reach. 
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